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ON THE MOTION OF TWO SPHEROIDS IN AN INFINITE LIQUID 
ALONG THEIR COMMON AXIS OF REVOLUTION. 

By Bibhutibhusan Datta. 

Introduction. 

Though the problem of the motion of two spheres in an infinite liquid 
along the line joining their centers has been completely solved by various 
investigators, the first writer to attempt the corresponding problem for 
two spheroids or ellipsoids is Prof. Karl Pearson.* His method does not, 
however, admit of further development and therefore, does not lead to the 
complete solution of the problem. In a previous paper, f I have shown how 
the problem can be completely solved in the case of two spheroids of revolu- 
tion of small ellipticities, the motion of the solids being along their common 
axis of revolution. The present paper deals with the more general case of 
the same problem in as much as no limitation has been imposed as regards 
the ellipticities. It will be seen that the success of the problem depends on 
certain transformation theorems for spheroidal harmonics which were not 
known before, though the corresponding theorems for spherical harmonics 
were given long ago by Bessel. 

All the results in this paper are believed to be new. 

§ 1. Lemma. 

1. Let ri, 0i, wi and r 2 , 2 , w 2 be two systems of spherical polar coordinates 
having their origins at 0i and 2 respectively and their polar axes in the 
directions of Oi0 2 and 2 0i; also let x u y h zi and x 2 , «/ 2 , s 2 be two parallel 
systems of rectangular cartesian coordinates with the same origins, the polar 
axes being taken for the axes of xi and xi\ let s denote the distance between 
0i and 2 . 

Since 



Pi = Vyl + z\ = pi 
we have 



* "On the Motion of Spherical and Ellipsoidal Bodies in Fluid Media," Part II (Quar- 
terly Journal of Pure and Applied Mathematics, Vol. XX). 
t Bulletin of the Calcutta Mathematical Society, Vol. VII. 
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when |s| > |xi + ipi cos 0\, we can expand the right-hand side of this 
equation in a series of spherical harmonics,* so that 

f n(CQS 6 2 ) (-1)" 1 1 f"[^ fe> , np f Sl + »Pi COS fl\ <*»&.«) 1 Jn 

where £ = s/ki and &i is any quantity. 

Following the notation of Sir W. D. Niven in his memoir " On Ellipsoidal 
Harmonics " (Phil. Trans., Vol. 128 (1891), A), let ©(2) denote the external 
harmonics of a system of confocal spheroids having the center at 2 and 
0\0 2 as the axis of revolution. Then ©(2) can be expressed in a series of 
spherical harmonics in the following way:| 



®w - (- w££yi.(±. 4;, £) { 



1+ " 



2(2n+3) 






2-4- (2ti + 3)(2n + 5) ' Hxl + yf+ z\' 

where D 2 stands for (a\ — c\)(d 2 jdxl), and a 2 , c 2 (a 2 > c 2 ) are the semi-axes 
of any one of the quadrics of the second family. 
Now 

/ JL A. iV-r-n" (2re ll IIn{x - 2 ' y2 ^ + 
Un \dx 2 ' dy 2 ' dz 2 )r 2 ~ { lj 2«7i! r^ 1 ' + 

and 

Hn(x 2 , 7/2, 2 2 ) = rlP n (C0S d 2 ) . 

Therefore 



„ (± ± AV = f _ n „(M!^ 

n \dx 2 ' dy 2 ' dz 2 )r 2 y ' 2 n n\ r n 2 



cos 2 ) 



d<? 



(X 2 + V2 COS l?) n+1 ' 

Hence if k\ denotes a\ — c\ we get 

* Vide Todhunter, "The Functions of Laplace, Lame and Bessel," p. 88. 

fVide Niven's memoir, p. 245; Hobson, "On the Evaluation of a Certain Surface 
Integral and its Application to the Expansion in Series of the Potential of Elipsoids," 
Proc. Lond. Math. Soc., Vol. XXIV, p. 91. 

t Niven's memoir, p. 236. 
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(2n + 1) 7T J L (*2 + *P2 cos #) re+1 

(»+!)(»+ 2) hi "I 

+ 2(2?i + 3) (x, + ip 2 cos t?)"+ 3 + " ' J m 

(2»)I 1 1 f" / ^+jp2_cos_g \ 7n 
~2»- 1 n!n!^ +1 xJ V " V k / 

Again 

® (2) = (2 W + i) »i fcF^l ro ?, (2m + 1) M — s; — ) d * 

R" (fa/tx)' ^ ...lorn 

A L *" 2(2w + 3) dt n+2 + J Wm w * 
Therefore we have 

(2n + 1) 1 1 1 r" ( x 2 + ip2Cosd 
2"n!n! AT^J, ^"V " *» 

= ^rw~J £ (2ra + 1)i H — fc, — )** 

V d n <j\ d n+2 "1 

X \_dT n + 2(2n + 3) dW 2 + " J Qm(<) ' (I) 
where <ri stands for fc 2 /&i. 

We shall now express this in terms of spheroidal coordinates. Supposing 
the spheroids are prolate, we write 

xi = ki cos 0i cosh rji = kifjL X \i, 

Pi = fci sin 6i sinh in = ft x (l - m 2 ) 1/2 (X 2 - 1) 1/2 ; 
then 

and 

^r g "( ^ 2+ t cos ' )^ =f "^" (X2) - 

Substituting in (I) we get 

PnfadQnfa) = E (2m + l)fii(m, n)P m 0*i)P m (\O, U) 

where 

Arfm.n) = (- l)"(2^ny! <r " +1 L^ + 2(2n+3)5^ + ■"JG-W- 

If however the spheroids are of the oblate or planetary form, a 2 — c\ is 
negative, so that we write y\ = c\ — a 2 , and the corresponding coordinates 
are 
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x\ — 71 cos 0i sinh t]i = 71M1X1, 
Pi = 71 sin 0! cosh m = 7l (l - M ?) 1/2 (X? + l) 1 ' 2 . 
Therefore (I) becomes 

P»G*«)Q»(iX«) = (- l) n+1 E (- l) m (2m+ l)fi;(m, n)P m {^)P m {ik^), (B) 

where 

N , 1X 2»re! ,n+irrf" , a', 2 d n+2 ~\ 

-,(«, ») = (- l)» (^+1)!-. L^' + 2(27+ 3) *P + J «»<« 

^i = 72/71 and h = 5/171. 

(A) and (5) are the lemmas that will be needed in the discussion of the 
problems that we shall presently take up. 

2. We can put the expressions for &i(m, re) and &i(m, re) into more 
compact forms as definite integrals. For 



d n Q m ( t) _ (- l)"re! f +1 P m (p)dp „ 



! f +1 Pm( 

Therefore 

2 n - 1 w!w! 
-,<», re) = (2 ^ ^ 



V r + T a " +1 1 (»+D(» +2) ^ -1 

X J- 1 Lr-^) ;r+ " i+ 2(2re+3) (T-^+i + ---J P ^)^ 

-•+1 



-sIMSt)*-^- 



Hence 

•+1 



*i(»», ») = \£* Qn (- - k * lP }p m (p)dp. 
Similarly we obtain 

*;o«, «) = IjT 1 Qn(^™)p m (p)d P . 



§ 2. Hydrodynamical Application. 

1. Let oi, Ci (ai > Ci) and a 2 , c 2 (a 2 > c 2 ) be the semi-axes of two ovary 
spheroids having their centers at 0\ and 2 respectively, O1O2 being the 
direction of motion. 

Let Mi be the velocity of 0i towards 2 , w 2 that of 2 towards d; also let 
s denote the distance between the centers 0i0 2 at any instant. Take Xi, jui, 
a and X 2 , ju 2 , « two systems of ovary spheroidal coordinates having their 
origins at 0\, 2 respectively, so that the surfaces of the given spheroids are 
Xi = X10 and X 2 = X 2 o. 

* Todhunter, loc. cit. 
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The problem before us is to find a function (j> satisfying the following 
conditions : 

(1) vV = 0, throughout the liquid, 

(2) cj> = 0, at infinity 

(3) rr- = — MifciPi(jui), when Xi = Xi , 

(Mi 

(4) rr- = - vtthPiQiz), when X 2 = X 20 , 

(7A2 

where fcf = a\ — c\ and h\ = a\ — c\. 

A value of <£ which satisfies the first two conditions is given by 

(5) <t> = E C^«P«(Mi)Q»(Xi) + ^„P„(M2)Q„(X 2 )]. 

11=1 

If by a proper choice of the arbitrary constants A n , B n , we can make it 
satisfy the other two conditions, then cj> will be the required solution. 
By the Lemma (A), near the surface of 0\, we have 

= Z [ A n PMQ n ^) + P n E (2m + l)S>i(m, n)P m ( Ml )P m (X 1 ) 1 

Substituting in the condition (3), we get 

- WiAaPiOl) = E ^«P«(mi)$>(Xio) 

»=1 L 

+ B E (2m + l)^i(m,n)P m ( M i)P;(Xi ) 1 , 

where Q'„(\io) and Pl(Xi ) stand for the values of (d/dXi)Qr>(Xi) and 
(d/dXi)P n (Xi) when Xio has been substituted for Xi. This must be true at 
every point on the surface of the spheroid. Hence we can equate to zero 
the coefficients of the various zonal harmonics of /ti. 
Equating the coefficients of Pi(jui) 

00 

- uxfcx = ^iQ'i(Xi ) + 3 E «i(l, n)B n . (6) 

«=i 

Equating the coefficients of P p (fj,i), p > 1, 

= A P Q' P Q, 10 ) + (2p + 1)P;(X 10 ) E &i(p, «)P„. (7) 

n-l 

We shall have two similar equations from the remaining surface condi- 
tion which can, however, be written down from symmetry. Thus 
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00 

- u 2 k 2 = BiQ[(\ 2 o) + 3 E «,(1, re)J„, (8) 

f or p > 1 

= 5 p Q;(X 20 ) + (2p + 1)P;(X 20 ) £ fi,(p, »)^„, (9) 

where &> 2 (m, re) can be written, by symmetry, from a>i(m, re), viz. 
w 2 (m, re) = ^ J i Q» ( — ^-- J P m (p)dp. 

Thus we have four equations (6) • • • (9) for the determination of the 
unknown constants A and B. 

2. The General Values of the Constants A and B. — Substituting the 
values of B\ and B n (n > 1) from the equations (8) and (9) in the equation 
(6), we get 

- mh = A 1 Qi(\ 10 ) - 3 ^Sr-v [ « 2 fc 2 + 3 E <S 2 (1, n)A n 1 

— 3 E (2m + l)wi(l, m) yfhr\ E -4rA(m, re), 

m=2 *£roVA20.J n-l 



or 



00 00 



E^4» E (2m + l)<5i(l, m)& 2 (m, re) n 7,.. 2 ° 

— iAiQi(Xio) = i«i&i — u 2 k 2 4 ' . -. 

Viv^2o; 

This equation can be written in the form 

00 

E OnmAn - i^iQi(Xio) = c, (m = 1, 2, 3, • • • ad inf) (1) 

»=i 

where 

0„ OT = (2m + l)«i(l, m)& 2 (m, re)P;(X 20 )/Ql(X 2 o), (2) 

and 

c = ^Wifci — m 2 ^2 S i(1, l)/Qi(X 2 o)- (3) 

In a similar manner, for the determination of B's, we get the equations 

E 0™5» - iJ5iQi(Xj ) = C, (m = 1, 2, 3, • • • ad inf.) (4) 

where 

C = (2m + 1)«S 2 (1, m)w!(m, n)P' m (K 10 )lQ' m (Kw), (5) 

and 

c' = Jw 2 fc 2 - wifci<M 2 (l, l)/Qi(Xio). (6) 

Thus we have finally an infinite number of equations (1) and (4) for 
the determination of an infinite number of unknowns, A as well as B. 
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The theory of the solution of the equations of this class has been worked 
out by Hill, Poincare, Koch, Toeplitz, Schmidt and others. Hence the 
values of the constants can be found so that the problems becomes deter- 
minate, the value of the potential function at any point being given by (5) 
Art. 1, §2. 

The general values of A as well as B, so determined will, however, be of 
little help for numerical calculations which will be necessary so as to get a 
physical idea of the state of liquid at any point. So we give a method of 
approximating to their values. 

3. Approximate Values of the Constants A and B. — We have 

2"n! V d n <s\ d n+2 1 

« l( m,n) = (- 1)» {2n + 1} j tf* |^ + 2{2n + 3) ^ 2 + ••■JQ„(0. 

Substituting the value of Q m (t) in series, viz. 

2™m!m! I" J_ (m +l)(ro+2) J_ 1 

Hm{t) (2m + 1) ! L t m+1 2(2m + 3) t m + 3 + " ' J 

and rearranging the terms we get 

_ 2 m+n m ! n ! (m + n) ! kTkp - 1 
Wl(m ' n) ~ (2m + 1) ! (2» + 1) ! ^ +1 

2 m + n mlnl(m+n + 2)1 fcrfc' 2 + 1 , 2 , 2 
+ (2m + 1) ! (2n + 1) ! 2(2w + 3) s m +"+ 3 ( x + 2J 

Thus the lowest order of «i(m, n) is 

/linear dimension \ m + n + l 
\ central distance / 

If the spheroids are so separated that we can neglect the terms of the order 

. ,. , I , we get from the equation (1), Art. 1, § 2 

central distance / 

A\= — n'r\ X \ > consequently Bi = — n , jX r, by symmetry. 
yi(Aio; V2IA20J 

We see generally that, for p > 1 

A » = " 2?+ 3!( ( 2^ 1)! ^(W ^ 5l + termS ° f higher ° rderS ' 
B » = ~ 2! ~lr^yr L! tS ^ ^ + terms of higher orders. 
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We proceed to approximate to the values of the A's and B's 
First Approximation. — 

2 Pl(Xx ) kV4 „ n 2 P' 2 (K 20 ) k\k\ 

M ~ 3 ^(x 10 j > ^ * 2 ~ " 3 q;(x 20 ") s~ M ' 

A = _ A P ^ ^ R R _ __ 1 P ^<>) #*I R 

15QJ(Xio) * 5 b 3 15Q:,(X 2 o)> *' 

Second Approximation. — 

Ml^l_ 2 1 _ k 2 kl 

1 ~ ~ QlCM + 3 Q[(Kio)Ql<M ~'s rU2 ' 
_ ji 2 k 2 2 1 _ kjztf 

1 ~ Qfcoj + 3 q;(x 2 o)q;(Xio) s 3 Ml> 

and so on. 

4. Verification of the Results Previously Obtained. — Neglecting terms of 

the order ( — - — . ,. ) and higher, we have found that 

\ central distance/ 

= " Q "^j Pi(Mi)Qi(Xi) - Q^j Pi(M»)Qi(X»). (1) 

Expressed in terms of spherical harmonics 

_ 2 J Pi(cos_0i) 3 P 3 (cos 00 



Pi(M2)Qi(X 2 ) = 3 «M "| H 4 &1 -|— + 

also 

2 f 2_ 34 

3 U 3 5 2 5 



«(»>--* 3+B+" 



Now if ei and e 2 be the eccentricities of the generating ellipses of the two 
spheroids respectively, we know 

C\ = 1/Xio, c 2 = I/X20, &i = aiCi, k 2 = a 2 e 2 . 

Substituting in the equation (1) and neglecting powers of c\ and e 2 higher 
than the second, we get as far as the required order 

.13/16 2N- P l( C0S ^l) 1 o. , 2 -P3(COS0l) 

, 1 s n fi 2n P i( cos ^) , 3 .-, 2 ^3(C0S 2 ) 

+ §M 2 Oj (1 - gCa) - ■ 2 + tV M 2«2 C 2 14 • 

r 2 r 2 



142 Datta: Motion of Two Spheroids. 

Now from the results of the previous paper, we have neglecting powers 
of ei and e 2 higher than the second 

*--**(&$)■ * — w(^)- 

Making allowance for the difference of notations between the two 
papers we have 

a = Oi(l - %el), b = o 2 (l - \e\), e t = \e\, e 2 = \e\. 

Therefore 

Ai \u x a\(l - |e?), Bi = - i«»aS(l - fe|). 

„. ,, _, , Al , /linear dimension V 

bince the terms or the order I — - — ; — 71— I are neglected, 

\ central distance/ ° 

^2 = 0, ^ 3 = - ^utalet, Ai = 0, 

JS 2 = 0, -B 3 = - -<hu 2 a\el, Bi = 0, • • • . 

Pi(cos'0O 



Further 



dsf 2 \ ri / 



3! : + 



Hence the equation (c) becomes up to the order of terms retained 

- 1 3/1 6 2 xPi(cos0i) 3 52 P 3 (cos0i) 
4> = iwiaf (1 - f e0 p 1- -fsuioiel -t — 

. 1 -vi 62 /i(cos« 2 ) 3 „ 2 P 3 (cos fl 2 ) 

'2 ~2 

so that the results of the two papers agree at least up to this order. Simi- 
larly it can be shown that the two results agree in their most general form. 

It may also be pointed out that, considered term for term, the series 
for the velocity potential function now obtained is more convergent than 
the series obtained in the previous paper. 

5. Conclusion. — In the preceding articles, I have studied the case when 
the two spheroids have the same common axis of revolution which coincides 
also with the direction of motion. In a similar way, can be investigated 
the case when the axes of revolution of the two spheroids are at right angles 
to the direction of motion and the transformation theorem for spheroidal 
harmonics that will be required can be obtained. 

University of Calcutta, 1919. 



